Observation of direct-photon collective flow in Au+Au collisions at y/s NN = 200 GeV 
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The second Fourier component V2 of the azimuthal anisotropy with respect to the reaction plane 
was measured for direct photons at midrapidity and transverse momentum (pr) of 1-13 GeV/c in 
Au+Au collisions at ^/s NN — 200 GeV. Previous measurements of this quantity for hadrons with 
Pt < 6 GeV/c indicate that the medium behaves like a nearly perfect fluid, while for pr > 6 GeV/c 
a reduced anisotropy is interpreted in terms of a path-length dependence for parton energy loss. 
In this measurement with the PHENIX detector at the Relativistic Heavy Ion Collider we find 
that for pr > 4 GeV /c the anisotropy for direct photons is consistent with zero, as expected if 
the dominant source of direct photons is initial hard scattering. However, in the pr < 4 GeV jc 
region dominated by thermal photons, we find a substantial direct photon V2 comparable to that 
of hadrons, whereas model calculations for thermal photons in this kinematic region significantly 
underpredict the observed V2- 

PACS numbers: 25.75.Dw 



Direct photons are produced in various processes dur- 
ing the entire space-time history of relativistic heavy ion 
collisions and, due to their small coupling, can leave the 
collision region without appreciable further interaction. 
This makes them a sensitive and direct probe of all stages 
of the collision, including initial hard scattering, forma- 
tion and evolution of the strongly interacting partonic 
medium, its transition to hadronic matter, and final de- 
coupling [l|, |2j. The transverse momentum (jpr) ranges 
populated by various production mechanisms overlap. 
However, azimuthal asymmetries tied to the event-by- 
event collision geometry provide useful additional infor- 
mation and a means to distinguish between sources of 
direct photons. In this paper we consider the second 
Fourier component (v2, often referred to as elliptic flow) 
of the event-by-event photon distribution in azimuth with 
respect to the reaction plane for minimum bias and se- 
lected centralities in Au+Au collisions. 

At higher pt (> 4 GeV/c) there are four fundamen- 
tal sources of direct photons, characterized by different 
V2 [H H[ • Photons from initial hard scattering (predom- 
inantly from qg — > 57 "gluon Compton scattering") are 
isotropic and so 1*2 = 0. Jet fragmentation photons have 
positive V2 since the energy loss of the originating par- 
ton is smaller in the reaction plane [4|. Jet-conversion 
photons where a hard scattered quark interacts with a 
thermal gluon in the medium and converts into a photon 
with almost equal px have negative vi 0|, because the 
average pathlength of the parton in the medium (pro- 
portional to the conversion probability) is larger out of 
the reaction plane than within. Finally, Bremsstrahlung 
photons are also emitted preferentially in the direction 
where the medium is thicker, leading to a negative V2 @|- 
Note that in this picture the azimuthal asymmetry of 
high pt photon production - while expressed in terms of 
V2~ reflects the pure geometry of the medium, not its dy- 



namics: it depends on the pathlength, not on the boost 
from the hydrodynamic pressure gradients. 

The picture is quite different in the low pr range 
(1 < pt < 4 GeV/c) dominated by thermal photons, as 
first measured in [5(, where bulk dynamics (expansion) 
plays an important role since it influences both the rate 
and azimuthal asymmetries of photon production 0, Q . 
It is now established that collectivity - which already 
exists in the partonic phase (strongly interacting Quark- 
Gluon Plasma, sQGP) - persists after transition into the 
hadronic phase and the resulting azimuthal asymmetries 
in particle production can be described by near-ideal hy- 
drodynamics. The expectation is that thermal radiation 
from both the sQGP and the hadronic phase will inherit 
the collective motion of the medium, i.e. will have a 
bona fide elliptic flow, positive V2 at low px [7] ■ The low 
Pt behavior of direct photon V2 puts constraints on the 
viscosity of the sQGP [6]. 

The PHENIX experiment has published the invariant 
yield as a function of pt for direct photons both via real 
photons and internal conversions of nearly real virtual 
photons [1, Q. In the 1 < pt < 4 GeV/c region, a 
substantial excess of direct photons was observed relative 
to scaling of p+p yields and has been interpreted in terms 
of thermal photon emission from the hot medium. An 
early attempt to infer V2 of direct photons from a ir° and 
inclusive photon V2 measurement performed in a limited 
Pt range has been published in [9(. In this Letter we 
present measurements by the PHENIX experiment [1(3] 
of V2 of 7r° and inclusive photons in a much extended 
transverse momentum (pt) range (up to 13 GeV/c) in 
y/s^" = 200 GeV Au+Au collisions. Also, at low pt 
the fraction i? 7 of direct over inclusive photons is now 
measured with much higher precision [j| than before Q , 
therefore, for the first time a meaningful extraction of 
the direct photon V2 itself is possible. 
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Data were taken in the 2007 run of the Relativis- 
tic Heavy Ion Collider at Brookhaven National Labora- 
tory. The analyzed sample includes ~ 3.0 x 10 9 mini- 
mum bias Au+Au collisions. Events were triggere d by 
the Beam-Beam Counters (BBC), as described in llj . 
which comprise two arrays of Cerenkov counters cover- 
ing 3.1 < \rj\ < 3.9 and 2-7T in azimuth in both beam 
directions (North and South). Event centrality was de- 
termined by the charge sum in the BBC. 

The event-by-event reaction plane (RP) has been de- 
termined by two detectors, the first being the BBC it- 
self. The RP resolution (effectively a dilution factor with 
which the observed «2 is normalized to obtain the true 
V 2 ) is defined as ct rp =< cos[2(* true - * RP )] > and it 
is established by comparing event-by-event the RPs ob- 
tained separately in the North and South detectors. The 
resolution is highest in the 20-30% centrality bin where it 
reaches a value of 0.4. For the 2007 data ta king period, 
a dedicated reaction plane detector (RXN) [12j was in- 
stalled covering 1.0 < < 2.8 and the full azimuth. The 
RXN is a highly segmented lead-scintillator sampling de- 
tector providing much better measurement (orp ^0.7) 
than the BBC, but it is closer to the central \rj\ < 0.35 
pseudorapidity region where V2 is measured, making it 
more sensitive to jet bias in those (rare) events where 
a high pt particle is observed. The 0.7/0.4 = 1.75 im- 
provement on the reaction plane resolution results is a 
1.75-fold improvement on point- by-point uncertainty. 

Inclusive photons were measured in the PHENIX elec- 
tromagnetic calorimeter 13J. Particles were identified 



(PID) and hadrons were rejected by a shower shape cut 
and a veto on charged particles using the Pad Cham- 
bers [14| . The remaining sample is collected for each pt 
range in histograms binned according to $ — <I/ RP where 
<f RP is the azimuth of the event-by-event reaction plane 
and established independently by the BBC and RXN. 
These distributions are then fit for each pr range with 
N a [1 + 2 v 2 cos{2($-* RP )}] to extract the raw ^' mcas 
coefficient for inclusive photons. As a cross-check of the 
fit value, another tijj' meas is also calculated from the av- 
erage cosine of the particles with respect to the reaction 
plane. While the PID eliminates virtually all hadrons 
above 6 GeV deposited energy (which might come from 
hadrons of any px above 6 GeV/c), a significant fraction 
of hadrons (up to 20% below 2 GeV deposited energy) 
survive the photon identification cuts. Since hadrons are 
known to have a large v 2 value, the observed «2 bs of in- 
clusive photons is obtained after correcting for hadrons 
as 



7,obs 



^^ynadr y/jymeas^ 



mcas\,,hadr 
v 1 



I _ _/yhadr jj^n 



where u hadr is the elliptic flow of hadrons and 
7V hadr /7V mcas is the fraction of hadrons in the sample 
surviving the PID cuts, as estimated from GEANT sim- 
ulations (20% at 2 GeV, 10% at 4 GeV and negligible 



above 6 GeV deposited energy). Finally the true i>2' mc 
for inclusive photons is obtained by dividing by the reac- 



tion plane resolution vj' 



7,obs 



/fRP- 



A large fraction of inclusive photons comes from 
hadron decays, predominantly from 7r° (~80%) and r\ 
(~15%), with a small fraction coming from p, u) and rj' 
decays, but only the ir° V2 is directly measured. The 
measurement of neutral pions and their v-i is described 
in detail in 0, EH- We assume that 7], u, etc. follow 
the same KEt scaling observed in hadrons [l6| where 
KEt — ti%t — m, Thus, v^^ipr) can be calculated for 
all hadrons from (pt)- For this we assume mr-scaling 
of hadron px spectra and establish a "hadron cocktail" 
using the measured yield ratios, similar to the one in [5j. 
This cocktail is the input of a Monte Carlo simulation 
to calculate the total v2' hs due to photons from hadron 
decays. The direct photon is then obtained using 

the R 7 (pt) "direct photon excess ratio" as 



7,dir 



i? 7 (p T W mC 



„7,bg 



where R^pt) = N inc (p T )/N b 3(p T ) with N inc = 7V moas - 
7V hadr , the number of inclusive photons, while N bg (pT) 
is the number of photons attributed to hadron decay. 
Values of R 7 (pr) above 5 GeV/c are taken from the real 
photon measurement with the PHENIX electromagnetic 
calorimeter Q, and below that from the more accurate, 
but pT-range limited internal conversion measurement of 
direct photons 0. 



TABLE I: Systematic uncertainties (Sx/x) contributing to the 
direct photon i>2 ,dlr measurement for minimum-bias collisions 
over two pr ranges. 



Contributing 
via 


Source 


PT 

1-3 GeV/c 


range 

10-16 GeV/c 


Type 




remaining hadrons 
t>2 extraction method 


2.2% 
0.4% 


N/A 
0.6% 


B 
B 




v 2 


particle ID 
normalization 
shower merging 


3.7% 
0.4% 
N/A 


6.0% 
7.2% 
4.0% 


B 
B 
B 


subtraction 


R ~, 


3.1% 


22% 


B 


common 


reaction plane 


6.3% 


6.3% 


C 



Sources of systematic uncertainties for representative 
Pt values are listed in Table U along with their characteri- 
zation: type A means point-by-point uncertainties which 
are uncorrelated with px, type B means uncertainties 
that are correlated (with pt) and type C is the over- 
all normalization uncertainty, moving all points by the 
same fraction up or down. Since the V2 measurement is a 
relative one (the azimuthal anisotropy is fit without the 
need to know the absolute normalization), the tt° and 
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FIG. 1: (Color online) (a,b,c) V2 in minimum bias collisions, 
using two different reaction plane detectors: (solid black cir- 
cles) BBC and (solid red squares) RXN for (a) 7r°, (b) inclu- 
sive photon, and (c) direct photon, (d) direct photon fraction 
R-y for (solid black circles) virtual photons -5] and (open blue 
squares) real photons [j| and (e) ratio of direct photon to tv° 
V2 for (solid black circles) BBC and (solid red squares) RXN. 
The vertical error bars on each data point indicate statistical 
uncertainties and shaded (gray and cyan) and hatched (red) 
areas around the data points indicate sizes of systematic un- 
certainties. 



inclusive photon v-i measurements are largely immune to 
energy scale uncertainties which are typically the domi- 
nant source of uncertainty in an absolute (invariant yield) 
measurement. The uncertainties on vi are dominated by 
the common uncertainty on determining grp and by un- 
certainties on particle identification. Uncertainties from 
absolute yields enter indirectly via the hadron cocktail 
(normalization) and more directly at higher px (where 
the real photon measurement is used) by the Ry(px) 
needed to establish the direct photon vi- Note that due 
to the way v\' Avc is calculated, once i? 7 is large, its rela- 
tive error contributes to the error on ' dlr less and less. 

Figure [T] shows steps of the analysis using the mini- 
mum bias sample, as well as the differences between re- 
sults obtained with BBC and RXN. The first v% of ir° and 
inclusive photons (v% ,v \ ' mc ) are measured, as described 



above (panels (a) and (b)). Then, using the v\" 9 of pho- 
tons from hadronic decays and the i? 7 direct photon ex- 



cess ratio, we derive the 



of direct photons (panel 



(c)). Panel (d) shows the Ry(px) values from the di- 
rect photon invariant yield measurements using internal 
conversion [fj and real Q photons, with their respective 

uncertainties. Panel (e) shows the ratio of v^ d%r jv^ . 
We observe substantial direct photon flow in the low px 



region (c) , commensurate with the hadron flow itself (e) . 
However, in contrast to hadrons, the direct photon V2 
rapidly decreases with px\ and starting with 5 GeV/c 
and above, it is consistent with zero (c). The rapid tran- 
sition from high direct photon flow at 3 GeV/c to zero 
flow at 5 GeV/c is also demonstrated on panel (e), since 
the 7T° v-i changes little in this region 
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FIG. 2: (Color online) (a,c,e) Centrality dependence of vi 
for (solid black circles) it , (solid red squares) inclusive pho- 
tons, and (b,d,f) (solid black circles) direct photons measured 
with the BBC detector for (a,b) minimum bias (c,d) 0-20% 
centrality, and (e,f) 20-40% centrality. For (b,d,f) the direct 
photon fraction is taken from 0] up to 4 GeV/c and from [sj 
for higher pr- The vertical error bars on each data point 
indicate statistical uncertainties and the shaded (gray) and 
hatched (red) areas around the data points indicate sizes of 
systematic uncertainties. 

A major issue in any azimuthal asymmetry measure- 
ment is the potential bias from where in pseudorapidity 
the (event-by-event) reaction plane is measured. At low 
Pt~ where multiplicities are high and particle production 
is dominated by the bulk with genuine hydrodynamic be- 
havior - there is no difference between the flow derived 
with BBC and RXN. However, at higher px we observe 
that the v-i values using BBC and RXN diverge, particu- 
larly for 7T° (panel (a) in Fig.QJ, less for inclusive photons. 
For direct photons (panel (c)) the two results are appar- 
ently consistent within their total errors, including the 
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error 8R 1 /R 1 (see Table H]) but it should be noted that 
Rj is a common correction factor in the v% measurements 
with both reaction plane detectors. 

Event substructure not related to bulk properties and 
expansion - most notably jets - can bias the reaction 
plane measurement, particularly at higher px and lower 
multiplicity. Observation of a high px particle practically 
guarantees the presence of a jet, which in turn modifies 
the event structure over a large rj range. The bias on the 
true event plane (with the bulk as its origin) is stronger if 
the overall multiplicity is small and if the rj gap between 
the central arm (where v% is measured) and the reaction 
plane detector is reduced. The bias in Fig. Q] is largest 
for 7T°, since high pt hadrons are always jet fragments. 
Inclusive photons are a mixture of hadron decay photons, 
inheriting the bias seen in ir° and the mostly unbiased 
direct photons, therefore, the difference between BBC 
and RXN is smaller. Finally, the bias is smallest (but 
nonzero) for direct photons, of which only a relatively 
small fraction (jet fragmentation photons) exhibit bias. 

Figure [2] shows v 2 for minimum bias and two central- 
ities as a function of transverse momentum for 7r°, in- 
clusive and direct photons. For reaction plane deter- 
mination the BBC is used because it is farthest from 
midrapidity where v 2 is measured. Despite the fact that 
there is a significant direct (thermal) photon yield at low 
Pt H j the tt° and inclusive photon v 2 is virtually identi- 
cal there. Note that the surprisingly large inclusive pho- 
ton v 2 is confirmed by the (so far preliminary) results 
with a completely different analysis technique [17j. For 
direct photons at low px we observe a pronounced pos- 
itive v 2 signal, increasing with decreasing centrality and 
comparable to the tt° flow, but then rapidly going toward 
zero at 5-6 GeV/c. Qualitatively this shape agrees with 
the prediction for very early thermalization times, 0.2- 
0.4 fm/c in [l8j], 0.2 fm/c and vanishing viscosity in Q, 
but both models severely underestimate the magnitude of 
the v 2 . The model in [19| combines somewhat later ther- 
malization time (0.6 fm/c) with partial chemical equi- 
librium in the hadronic phase, reproducing the shape, 
but missing the magnitude of the observed v 2 at low px ■ 
While such large direct photon v 2 in principle could be 
attributed to a dominant production mechanism at the 
later stage when bulk flow is already developed, simul- 
taneously explaining the large values of v 2 at ^2 GeV/c 
and its vanishing above 5 GeV/c remains a challenge to 
current theories. 

Figure |3] shows the high px integrated v 2 (px > 
6 GeV/c) for tt° and photons (inclusive and direct) as a 
function of centrality. The low N par t behavior is strongly 
influenced by the location in pseudorapidity of the re- 
action plane detector. The tt° v 2 is comparable to other 
hadrons and is higher than the inclusive photon v 2 , which 
is diluted by direct photons. The two direct photon v 2 
measurements (panel (c)) are consistent with zero (and 
each other) at all centralities within their total system- 
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FIG. 3: HighpT (pt > 6 GeV/c) integrated V2 vs N pal t for (a) 
7T°, (b) inclusive photon, and (c) direct photon. Results are 
shown with both reaction plane detectors: (solid black circles) 
BBC and (solid red squares) RXN. Each point represents a 
10% wide centrality bin from 60-0%. The vertical error bars 
on each data point indicate statistical uncertainties and the 
shaded (gray) and hatched (red) areas around the data points 
indicate sizes of systematic uncertainties. 

atic errors. While zero v 2 would be expected if initial 
hard scattering is the dominant (sole considered) source 
of photons, it should be pointed out that the typical con- 
tribution from jet-conversion only would be v 2 ~ —0.02 
and from fragmentation v 2 < 0.01 weighted with the frac- 
tion of photons coming from these specific processes [3, 0] ■ 
Currently the experiment is not sensitive to their nega- 
tive/positive contributions to v 2 . 

In conclusion, PHENIX has measured v 2 of 7r°, inclu- 
sive and direct photons in the 1 < px < 13 GeV/c range 
for minimum bias and selected centralities in y?s NN = 
200 GeV Au+Au collisions. At higher p T (> 6 GeV/c) 
the direct photon v 2 is consistent with zero at all cen- 
tralities, as expected if the dominant source of photon 
production is initial hard scattering. However, the ex- 
perimental uncertainties are currently about a factor of 
2 higher than the predicted (small) positive and negative 
contributions from fragmentation and jet conversion pho- 
tons, respectively. In the thermal region (px < 4 GeV/c), 
a positive direct photon v 2 is observed which is compara- 
ble in magnitude to the ir° v 2 and consistent with early 
thermalization times and low viscosity, but its magnitude 
is much larger than current theories predict. 
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